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Some aspects of heavy flavor probes of quark-gluon plasma (QGP) including quarkonia and 
open heavy flavor D- and B-mesons, aimed to study the properties of QGP, are discussed in 
this lecture note. 

1 Introduction 

Lattice QCD (LQCD) calculations predict that at a critical temperature Tc ~ 170 MeV, cor- 
responding to an energy density Ec ~ 1 GeV/fm^, nuclear matter undergoes a phase transition 
to a deconfined state of quarks and gluons, called Quark-Gluon Plasma (QGP). At the modern 
collider facilities such as CERN supersynchrotron (SPS) (the nucleon-nucleon (NN) centre-of- 
mass energy for collisions of the heaviest ions is ^/s = 17.3 GeV), Relativistic Heavy Ion Col- 
lider (RHIC) at Brookhaven (^ = 200 GeV), and Large Hadron Collider (LHC) at CERN 
{^/s = 5.5 TeV), whose heavy-ion program will start soon, heavy-ion collisions are used to attain 
the energy density, exceeding Sc- This makes the QCD phase transition potentially realizable 
within the reach of the laboratory experiments. The objective is then to identify and to assess 
suitable QGP signatures, allowin g to study the properties of QGP. To that end, a variety of 
observables (probes) can be used^^^. Further we will be mainly interested in heavy-flavour 
probes of QGP, i.e., utilizing the particles having c- and b-quarks. 

A special role of heavy Q = (c, b) quarks as probes of the medium created in heavy-ion 
collision (HIC) resides on the fact that their masses {rric ~ 1.3 GeV, ~ 4.2 GeV) are signifi- 
cantly larger than the typically attained ambient temperatures or other nonperturbative scales, 
niQ ^ Tc, A.QCD = 0.2 GeV^. This has several implications: (i) The production of heavy quarks 
is essentially constrained to the early, primordial stages of HIC. Hence, heavy quarks can probe 
the properties of the dense matter produced early in the collision, (ii) Thermalization of heavy 
quarks is "delayed" relative to light quarks. One could expect that heavy quarks could "thermal- 



ize" to a certain extent, but not fully on a timescale of the lifetime of the QGP. Therefore, their 
spectra could be significantly modified, but still would retain memory about their interaction 
history, and, hence, represent a valuable probe, (iii) RHIC, and especially LHC experiments 
allow to reach very low parton momentum fractions x, where gluon saturation effects become 
important. Heavy quarks are useful tools to study gluon saturation, since, due to their large 
masses, charm and bottom cross sections are calculable via perturbative QCD and their yield is 
sensitive to the initial gluon density. 

The heavy-flavor hadrons we will be interested in include: 1) open charm D = (cq) and 
open beauty B = {hq) mesons composed of a heavy quark Q = (c, 6) and a light antiquark 
q = {u,d). These mesons could be sensitive to the energy density of the medium through the 
mechanism of in-medium energy loss; 2) hidden charm [charmonia= (cc)] and hidden beauty 
[bottomonia=(65)] mesons (called collectively heavy quarkonia) being the bound states of the 
charm quark-antiquark, or bottom quark-antiquark pairs, respectively. Heavy quarkonia could 
be sensitive to the initial temperature of the system through the dissociation due to color 
screening of the color charge that will be discussed later. 

For detecting heavy flavor hadrons, different decay channels are used. At LHC, the ALICE 
experiment for detecting open charm will use the reconstructed hadronic decays like — )■ 
K~Tr^, —7- i^^vr^vr"'", etc., for detecting open beauty ALICE will use inclusive single lepton 
decays B ^ e + X (midrapidity) , B ^ fi + X (forward rapidity) and inclusive displaced 
charmonia decays B J/ipi-^ e+e~) + X. At RHIC, in PHENIX and STAR experiments 
the measurement of the spectra of open heavy flavors is based on the measurement of the 
spectra of heavy flavor (HE) electrons and positrons [(e"'' + e~)/2] from the semileptonic decays 
like — >• K~e'^i'e, — )• K^e'^i'e, etc. These measurements are based on the fact that 
the decay kinematics of HE electrons/positrons largely conserves the spectral properties of the 
parent particles. Besides, the STAR experiment has the capability to directly reconstruct open 
heavy flavor mesons through the hadronic decay channels {D'^ — )■ Ktt, — )• i^vrvr, etc.). At 
low pt, the STAR experiment also uses heavy flavor decay muons to provide open heavy flavor 
measurements. In all of the above experiments, quarkonia are detected through their dilepton 
decays QQ — )• e+e~ (midrapidity), QQ — )■ n'^n" (forward rapidity). 



2 Heavy Flavor Probes of QGP: Quarkonia 

We begin the discussion of heavy flavor probes of QGP with heavy quarkonia. The question we 
would like to address is: What could happen with quarkonium yields in HIC if QGP is really 
formed? In the above mensioned experiments, quarkonia are detected in the dielectron channel 
at midrapidity and in the dimuon channel at forward rapidity. Eor example, in the dimuon 
channel, signals of about 7 • 10^ J ftps and 10^ Ts are expected in the most central Pb-Pb 
collisions at LHC / ALICE for one year of data taking at nominal luminosity (see Fig.[I]from'21). 
If quarkonia will be completely suppressed by QGP, these signals will disappear, otherwise, in 
the case of enhancement of quarkonium yields, the peaks will be more pronounced as compared 
to the signals obtained in pp collisions after scaling them with the corresponding number of 
independent binary NN collisions. 

Let us note that if the deconfinement phase transition really takes place in a dense medium 
created in HIC, then a color charge in the QGP will be screened analogously to the Debye 
screening of an electric charge in the electromagnetic plasma. As a result of color Debye screening 
of the heavy quark interaction in QGP, the binding energy of a bound state decreases and one 
can expect that this would lead to the suppression of quarkonium yields in HIC. This idea was 
first suggested by H. Matsui and H. Satzl^who predicted that color Debye screening will result 
in the suppression of J/ip meson (cc in the ^5*1 state, M = 3.097 GeV) yields. After that, the 
J/ijj suppression was considered as one of the key probes for the QGP formation in heavy ion 
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Figure 1: Left: Dimuon invariant mass spectra (after subtraction of the uncorrelated background) in central 
Pb-Pb collisions at LHC/ALICE for one year of data taking (10^ s) and luminosity 5 • lO^^cm^^s^^ in the J/ijj 
mass region (PYTHIA calculations). Right: Same as on the left, but in the T mass region. 



collisions. J/^|J is especially promising because of the large production cross-section and dilepton 
decay channels which make it easily detectable. 

However, soon it was realized that besides melting J/^lJ mesons in QGP due to the screening 
of the color charge, there are also a few competing mechanisms which could explain the suppres- 
sion of J/^l) production in heavy- ion collisions. These mechanisms are referred to as cold nuclear 
matter (CNM) effects. The first CNM effect is the absorption of J/?/^ by nuclear fragments from 
colliding nuclei. Let us consider, e.g., the proton-nucleus collision. Once produced in the hard 
primary parton processes, J ftp has to cross the length L of nuclear matter, before exiting the 
nucleus, and, when traversing nuclear matter, it can be absorbed by forthcoming nucleons of a 
nucleus. The production cross section of J/V' is p-A collision can be parameterized as 

-ii' = A4l^e-<i^^\ (1) 

where Upp^ is the production cross section of J/ij) in pp collisions and <7^/^ is the nuclear 
absorption cross section. From the global fit to the data on charmonium production in p-A 
collisions the value of a'^J^^ can be extracted, in particular, at SPS (NA50 experiment) it was 

obtained that aH'^' = 4.2 ± 0.5 mbEl. 

The second CNM effect is related to shadowing of low momentum partons. This means the 
depletion of low momentum partons in nucleons bound in nuclei as compared to free nucleons. 
This effect can be accounted for in terms of the modification of the parton distribution functions 
in nucleon within the nucleus with respect to the parton distribution functions in a free nucleon: 



Rt{x,Q^)= %y:^J, <1, i = q.,qsea,g. (2) 



Here "i" denotes valence quarks, sea quarks, and gluons, x is the parton momentum fraction, 
is the momentum transfer squared. At high energies, J/^/^s are dominantly produced through 
the gluon fusion, and the J/ip yield is therefore sensitive to gluon shadowing. The underlying 
idea explaining the occurrence of gluon shadowing is that the gluon density strongly rises at 
small X to the point where gluon fusion, gg — t- g^ becomes significant. In the case of proton- 
nucleus and nucleus-nucleus collisions, where nuclei with large mass number A are involved, the 
nonlinear effects are enhanced by the larger density of gluons per unit transverse area of the 
colliding nuclei. A direct consequence of nuclear shadowing is the reduction of hard-scattering 
cross sections in the phase-space region characterized by small-x incoming partons. For gluons, 




Figure 2: J/xp yields normalized by Drell-Yan process, as a function of the nuclear thickness L, as measured at 

the SPS. 

e.g., shadowing becomes important at x < 5 x 10"^, and, hence, is relevant for the conditions of 
RHIC and LHC. Note however that the strength of the reduction is constrained by the current 
experimental data only for x > 10~^. 

After discussing CNM effects on the J/-^ production, let us consider the J/il) production at 
CERN SPS experiments. Fig.[2]shows the J/ijj production cross section, where J /ip are registered 
through their decay in the dimuon channel, normalized to the cross section of the Drell-Yan (DY) 
proces^. The nuclear thickness parameter L parameterizes the number of participating nucleons 
and increases with centrality. At SPS energies, DY process is not modified by the medium and 
its production cross section scales linearly with the number of binary NN collisions. Hence, this 
process can be used as unsuppressed reference process for the J/ip production. The current 
interpretation of the obtained results is the following. The effects of cold nuclear matter are 
quantified in proton-nucleus collisions, where the energy density is not enough to reach the 
critical value for the formation of QGP. There is no anomalous suppression, beyond CNM 
effects, of the J /ip yield in the collisions of light ions of sulfur [A = 32) with the heavy ions 
of uranium [A = 238). However, there is the anomalous suppression beyond the CNM effects 
in the central collisions of lead nuclei when the number of participating nucleons is the largest, 
and, to a certain extent, in the central collisions of indium nuclei {A = 115)^. This could be 
considered as the signature of melting the charmonium states in quark-gluon plasma. At this 
point, it is also worth to note that some level of the J /ip suppression could originate with the 
reduced feed-down to J/ip from excited charmonium states (V^'iXc)) which melt just above the 
QGP transition temperature. One needs further to investigate what is really suppressed, the 
directly produced J/V's or the J/ips coming from the heavier charmonium states. 

Let us consider the J/ij: production at RHIC experiments, when the nucleon-nucleon center- 
of-mass energy for collisions of the heaviest ions increases by a factor 10 compared to the CERN 
SPS experiments. The J/tp suppression can be characterized by a ratio called the nuclear 
modification factor 

(fNjU/dpTdy 

obtained by normalizing the J/ ip yield in nucleus-nucleus collision by the J/ ip yield m. p + p 
collision at the same energy per nucleon pair times the average number of binary inelastic NN 
collisions. This ratio characterizes the impact of the medium on the particle spectrum. If heavy 
ion collision is a superposition of independent N^-ou inelastic NN collisions, then Rab = 1, 
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Figure 3: Left: J/il)'s Raa for Au-Au collisions at midrapidity compared to a band of theoretical curves for the 
breakup values found to be consistent with the d-Au data. Both EKS and NDSG shadowing schemes are included. 

Right: Same as on the left but at forward rapidity. 



whereas Rab < 1 (Rab > 1) corresponds to the case of the J/^p suppression (enhancement). As 
we discussed aheady, at first, it is necessary to clarify the role of CNM effects on J/ip production. 
At RHIC, CNM effects are studied in collisions of light deuteron and heavy gold nuclei, when 
the energy density reached in the collision is not enough for the formation of QGP. At RHIC 
energies, shadowing of partons is important and in the model calculations is implemented in two 
shadowing schemes for the nuclear parton distribution functions, the EKS mo deM and NDSG 
mo The J/ip break-up cross sections obtained for two shadowing schemes from the best fit 
to data are ^breakup = 2.8l|f mb (EKS) and abreakup = 2.611^ mb (NDSG)[I2] (in Ref.M the 
term " break-up cross section" is used instead of the term " absorption cross section" ) . Although 
these values are consistent, within large uncertainties, with the corresponding value obtained at 
CERN SPS, a recent analysisEl, hows that, in fact, the level of J/^ CNM break -up significantly 
decreases with the collision energy. 

Let us now consider charmonium production in heavy ion collisions at RHIC. Fig.[3]shows the 
PT^-integrated J/ip nuclear modification factor obtained in Au-Au collisions at RHIC/PHENIX 
experiment as a function of centrality, parametrized by the number of participating nucleons 
at mid- (left panel) and forward (right panel) rapidity ^1^. The PHENIX data are shown by 
blue symbols. The Raa approaches unity for the peripheral collisions (small Npart) and goes 
down to approximately 0.2 at most central collisions (large Npart) ■ To see the level of the 
anomalous suppression beyond the cold nuclear matter effects it is necessary to extrapolate the 
CNM effects obtained in d-Au collisions to the Au-Au collisions within the given shadowing 
scheme and the J/ip break-up cross section. The results are shown by black and red curves with 
the corresponding error bands. It is seen that J /ip production is significantly suppressed beyond 
CNM effects at forward rapidity (right) and suppression is less pronounced at midrapidity (left) 
in most central Au-Au collisions. 

However, not all is still clear. Measurements of the J/^ suppression by PHENIX collabo- 
ration at RHIC lead to some surprising features. Fig. [J] shows compiled data for the nuclear 
modification factor obtained in CERN SPS and RHIC PHENIX experiments. There are two 
surprising results in these measurements. First, the mid rapidity suppression in PHENIX (the 
red boxes) is lower than the forward rapidity suppression (blue boxes) despite the experimental 
evidence that the energy density is higher at midrapidity than at forward rapidity, and, hence, 
one could expect that at midrapidity the J /tps will be more suppressed due to higher density of 
color charges. Secondly, the nuclear modification factor Raa at midrapidity in PHENIX (red 
boxes) and SPS (black crosses) are in agreement within error bars, a surprising result considering 
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Figure 4: J/ip nuclear modification factor for the most energetic SPS (Pb-Pb) and RHIC (Au-Au) collisions, as 

a function of the number of participants Npart- 



that the energy density reached at RHIC is larger than the one reached at SPS. This indicates 
that at RHIC energies additional mechanisms countering the suppression, could be operative. 

Let us consider possible explanations of the above features. 1. Regeneration of J /ips in the 
hot partonic phase from initially uncorrelated c and c quarks (quark coalescence model). If to 
compare the J/ip suppression pattern at RHIC and SPS, J/ips could be indeed more suppressed 
at RHIC than at SPS, but then regenerated during (or at the boundary of) the hot partonic 
phase from initially uncorrelated c and c quarks. If to compare the results at RHIC (midrapidity 
vs. forward rapidity), at midrapidity, due to higher energy density, there are more c and c to 
regenerate than at forward rapidity that could explain the stronger suppression at forward 
rapidity. Note that the total number of initial cc pairs is larger than 10 in the most central 
Au-Au collisions. Certainly, if regeneration is important at the RHIC conditions, it will be even 
more important at the LHC conditions where more than 100 cc pairs is expected to be produced 
in the central Pb-Pb collisions. 2. J /ip production could be more suppressed at forward rapidity 
due to the nuclear shadowing effects. Standard gluon shadowing parametrizations do not tend 
to produce such an effect but they are poorly constrained by the data and further saturation 
effects are not excluded^. 

To show the complexity of the problem, let us consider some theoretical models for the char- 
monium production, which quite satisfactory describe the RHIC data but whose predictions 
for LHC are drastically different. First, in the statistical hadronization model (SHM)^^!, it is 
assumed that: 1. All heavy quarks (charm and bottom) are produced in primary hard colli- 
sions and their total number stays constant until hadronization. 2. Heavy quarks reach thermal 
equilibrium in the QGP before the chemical freeze-out (hadronization). 3. All quarkonia are pro- 
duced (nonperturbatively) through the statistical coalescence of heavy quarks at hadronization. 
Multiplicities of various hadrons are calculated with the grand canonical ensemble. The gener- 
ation of J/ip proceeds effectively if c, c quarks are free to travel over large distances implying 
deconfinement . 

Fig. [5] shows the rapidity dependence of the nuclear modification factor, obtained in this 
model and the comparison with the rapidity dependence at PHENIX for two centrality bins. 
Two theoretical curves correspond to two fitting procedures, with one and two Gaussians of 
the J/^ data in pp collisions. In both cases, calculations reproduce rather well (considering 
the systematic errors) the Raa data. The model describes the larger suppression away from 
midrapidity. The maximum of Raa at midrapidity in this model is due to the enhanced gener- 




Figure 5: Rapidity dependence of Rjl^ for two centrality classes in the statistical hadronization model. The data 
from the PHENIX experiment (symbols with errors) are compared to calculations (lines, see text). 




Figure 6: Centrality dependence of the R/^ at midrapidity, according to the statistical hadronization model. 



ation of charmonium around midrapidity, determined by the rapidity dependence of the charm 
production cross section. The centrahty dependence of Raa at y = is shown in Fig. [6l The 
model reproduces quite well the decreasing trend with centrality seen in the RHIC data. Fig. [6] 
also shows the prediction of the model for the LHC. At much higher LHC energies, the charm 
production cross section is expected to be larger by about an order of magnitude. As a result, 
a totally opposite trend as a function of centrality is predicted, with Raa exceeding unity for 
central collisions. 

Let us consider the comovers interaction model (CIM)^!^. This model does not assume 
the deconfinement phase transition. Anomalous suppression of J/ip (beyond CNM effects) is 
the result of the final state interaction of the cc pair with the dense medium produced in the 
collision (comovers interaction). The model consistently treats the initial and final state effects. 
The initial state effects include: 1) nuclear absorption of the pre-resonant cc pairs by nucleons 
of the colliding nuclei, 2) consistent treatment of nuclear shadowing for hard production of 
charmonium. The final state effects include absorption of the cc pairs by the dense medium 
created in the collision (interaction with comoving hadrons produced in the collision). The 
model does not assume thermodynamic equilibrium and, thus, does not use thermodynamic 
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Figure 7: Results for J /il> suppression in Au-Au collisions at RHIC at mid- (left panel), and at forward rapidities 
(right panel) in the comovers interaction model. The solid curves are the final results. The dashed-dotted ones 
are the results without recombination (C = 0). The dashed line is the total initial-state effect. The dotted line 
in the right panel is the result of shadowing. In the left panel the last two lines coincide. 



concepts. The density of charmonium is governed by the differential rate equation 



dn 



dr 



-crco[nco{b, s, y)nj/^,{b, s, y) - ndb, s, y)nc{b, s, y)], 



(4) 



supposing a pure longitudinal expansion of the system and boost invariance. In Eq. ([4]), rico is the 
density of comovers, which is found in the dual parton mo 

delttSl 

together with the proper shad- 
owing correction, aco is the cross section of J/tl^ dissociation due to interactions with comovers, 
taken such as to reproduce the low energy SPS experimental data (with aco = 0.65 mb). The 
first term on the right describes dissociation of charmonium due to interaction with comovers. 
The second term describes the recombination of charmonium and is proportional to the product 
of densities of charm quarks and antiquarks. The important feature of the CIM is that recom- 
bination of c-c quarks proceeds only locally, when the densities of quarks and antiquarks are 
taken at the same transverse coordinate s. This is different from the recombination in the SHM, 
where recombining quarks can be separated by large distance that implies deconfinement. The 
effective recombination cross section in the CIM is equal to the dissociation cross section due to 
the detailed balance. The results for the centrality dependence of the Raa are shown in Fig. [71 
At midrapidity, the experimental data are well reproduced by full theoretical calculations (solid 
curve) taking into account nuclear shadowing, dissociation by comovers and recombination from 
charm quark and antiquark pairs. At forward rapidity, the results also well agree with data, in 
particular, the J/ip suppression at forward rapidity is somewhat larger that the suppression at 
midrapidity. Fig. [8] shows the predictions of the model for LHC. The parameter C encodes the 
recombination from c-c pairs and is vanishing in the absence of recombination. Although the 
density of charm grows substantially from RHIC to LHC, the combined effect of initial-state 
shadowing, absorption and comovers dissociation overcomes the effect of parton recombination. 
This is in sharp contrast with the predictions of the statistical hadronization model where a 
strong enhancement of the J/'i/' yield with increasing centrality was predicted. 

Thus, the J/tl^ suppression is an important characteristic to search for QGP. But if J/ip 
anomalous suppression (beyond CNM effects) was observed in HIC, there are a few competing 
mechanisms to explain that: 1. Charmonium is dissociated due to the genuine color screening 
in the deconfined medium. 2. Charmonium is dissociated through interactions with comoving 
hadrons in the medium formed in HIC. How to differ these mechanisms? In comovers inter- 
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Figure 8: Results for J/tl) suppression in Pb+Pb at LHC (yi = 5.5 TeV) at midrapidity for different values of 
the parameter C in the comovers interaction model. The upper line is the suppression due to initial-state effects 

(shadowing) . 

action model, the anomalous suppression sets in smoothly from peripheral to central collisions 
rather than in a sudden way when the corresponding dissociation temperature in the deconfined 
medium is reached. At SPS, current experimental errors still do not allow to disentangle these 
two mechanisms. However, even if the color screening mechanism is dominating, it is unclear 
what is really melted, directly produced J/^l's or originating from the feed-down of less bound 
charmonium states, Xc J/i^ + i^' J/i^ + which have lower dissociation tem- 
peratures. At RHIC and LHC conditions, the feed-down from B-meson decays becomes also 
important. Recombination enhances J/ip production and much complicates the picture but its 
effect may be different depending on whether the deconfinement phase transition happened or 
not. It is even possible that after all J ftps were melted in QGP they can be statistically regen- 
erated at hadronization. True operative mechanisms of J/^lJ production can be established only 
after studying all important dependences (from centrality, rapidity, collision energy ^/s,...) of 
all relevant observables with sufficient accuracy. 

So far we considered the suppression (or enhancement) patterns for the J/^lJ production. 
Now let us briefly consider the perspective for bottomonia. One could expect that bottomonium 
production might be easier to understand than charmonium production due to the following 
reasons. 1. Since less than one bb pair is produced in one central Au-Au collision, the regen- 
eration is negligible at the conditions of RHIC. Besides, only about 5 bb pairs are expected to 
be produced in a single central Pb-Pb collision at LHC. Hence, regeneration should play much 
less role in the beauty sector than in the charm sector. 2. Having higher masses, bottomonia 
originate from higher momentum partons and will less suffer from shadowing effects. 3. The 
absorption cross section for T is by 40 — 50 % smaller than the corresponding cross section for 
J/ijj and ip'. 

These features should ease the separation of the anomalous suppression in the T's family. 
3 Heavy Flavor Probes of QGP: Open Heavy Flavor Mesons 

What qualitative effects could one expect to obtain when probing the dense matter by heavy 
quarks (charm or bottom)? As well known from electrodynamics, the bremsstrahlung off an 
accelerated heavy quark Q is suppressed by the large power of its mass ~ (mg/mg)^ as compared 
to light quarks. Therefore, gluon radiation off heavy quarks (i.e., radiative energy loss) is much 
suppressed relative to light quarks. Consequently, one could expect a decrease of high pT 
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Figure 9: Nuclear modification factor (upper panel, central Au-Au) and elliptic flow (lower panel, minimum- 
bias Au-Au) of non-photonic electrons at RHIC, compared to theory. The band corresponds to the Langevin 
simulations based on an expanding fireball with effective heavy quark resonance interactions [22]. 



suppression and of the elliptic flow coefficient f from light to charm to bottom quarks. Or, 
that the energy loss and coupling to matter of heavy quarks is smaller than for light quarks as well 
as that the thermalization time for heavy quarks is longer than for light quarks. Due to the above 
features, one should observe a pattern of gradually increasing Raa when going from the mostly 
gluon-originated light-flavor hadrons {h^ and tt^) to D to B mesons: R^a ^ ^AA ^ -^AA^^^ 
(gluons lose more energy than quarks since gluons have a higher color charge) . The enhancement 
above the unity of the heavy-to-light ratio R^^J^ = R^j^/R\^ probes the color charge dependence 

of the parton energy loss while the ratio R^^^ = RaaZ-^aa Probes the mass dependence of the 
par ton energy loss. 

Let us now consider what the experiment tells us about the open heavy flavor px suppres- 
sion and elliptic flow. As mentioned earlier, at RHIC, open heavy flavors are studied through 
measurements with the electrons and positrons originating from the semileptonic decays of D- 
and i?-mesons. Fig. [9] shows the nuclear modification factor and elliptic flow coefficient for 
HF ele ctrons as functions of px, obtained in central collisions of gold nuclei at RHIC (closed 
circles)!!^. In contrast to the above expectations, the results for R^^ show a strong suppres- 
sion of HF decay electrons at px > 2 GeV/c, approaching at high px the level of suppression 
for tt". This evidences that produced medium is quite dense for heavy quarks to lose energy 
as efficiently as light quarks do. The measurement of elliptic flow gives rather large value for 
V2^. This means that HF electrons are involved in a collective motion being indicative of the 
collective flow of their parent particles as well. 

In Fig. [9l the results of some model considerations are also shown. The best description is 
provided by the model assuming the Brownian motion of heavy quarks within the framework 
of Langevin dynamics ^^Sl. Let us consider the basic assumptions of this theory. Firstly, the 
thermal heavy quark momentum p'^ ~ mT (in nonrelativistic approximation) is much larger 
than the typical momentum transfer ~ from a thermal medium to a heavy quark. Hence, 
the motion of a heavy quark in the QGP can be represented as the Brownian motion, which 
can be described using the Langevin equation. Secondly, heavy quark loses its energy in elastic 
scattering processes with light partons. Besides, as evidenced from calculations of heavy and 
light meson correlators within LQCD, in QGP the D- and 5-meson like resonant states exist 
up to T < 2Tc. Rescattering on these resonant states plays an important role in thermalizing 
heavy quarks. Thirdly, in order to get the spectrum of HF electrons, c- and 5-quarks are to be 



hadronized into D- and i?-mesons via quark coalescence (at low pt) and fragmentation (at high 
Pt)- 

The analysis shows that resonance scattering decreases nuclear modification factor and 
increases azimuthal asymmetry Heavy- light quark coalescence in subsequent hadronization 
significantly amplifies and increases R^^, especially in the pT — 2 GeV/c region. The 
contribution from B-mesons to R^^ and is estimated providing full calculations with c+b 
quarks and with only c quarks. The result is that the B-meson contribution increases R^^ and 
decreases V2^, and becomes important above — 3 GeV/c. Thus, one can conclude, that the 
combined effects of coalescence of heavy quarks Q with light quarks q, and of the resonant heavy- 
quark interactions are essential in generating strong elliptic fiow of up to 10%, together 
with strong suppression of heavy fiavor electrons with R^^ about 0.5. 
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